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ABSTRACT 


Damping  of  welded  structures  is  a  subject  of  great  interest  and  application  for  the 
navy  as  relates  to  ship  shock  survivability  and  acoustic  transmission  of  ship  noise.  The 
purpose  of  this  research  is  to  study  the  effects  of  welding  on  damping.  A  generic  model 
of  a  warship’s  hull  structure  was  used  to  study  damping  effects.  The  model’s  natural 
frequencies  and  mode  shapes  were  calculated  using  a  finite  element  model  prior  to  model 
testing.  The  frequency  response  and  natural  frequencies  of  the  model  were  determined 
experimentally  by  exciting  the  model  and  measuring  the  response  throughout  the 
structure  using  Frequency  Response  Functions  (FRF’s).  The  results  were  compared  with 
the  finite  element  modeling.  The  damping  ratio  of  the  model  in  relation  to  position  from 
excitation  was  calculated  using  the  half-power  point  method  and  then  a  more  detailed 
analysis  of  frequency  dependent  damping  versus  position  was  made  using  modal 
parameter  extraction  using  the  Complex  Exponential  Method. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Damping  of  ship  structures  is  a  subject  of  great  interest  and  application  for  the 
navy  as  it  relates  to  ship  shock  survivability  and  acoustic  transmission  of  ship  noise.  The 
concept  of  damping  can  be  thought  of  as  the  culmination  of  various  energy  dissipation 
mechanisms  that  remove  mechanical  energy  from  a  vibrating  system.  The  energy 
dissipation  mechanisms  can  be  broadly  categorized  into  two  categories,  hydrodynamic 
effects  and  structural  damping  of  the  ship  hull.  While  some  study  has  gone  into  the 
hydrodynamic  effected  hull  damping,  little  study  has  gone  into  structural  damping  in  the 
ship  hull.  The  ship  system  has  many  proposed  energy  dissipation  sources  such  as  welded 
joints  and  stiffeners,  long  cable  trays,  hangers,  snubbers,  etc.  This  thesis  focuses  on  the 
effect  of  welds  in  the  ship’s  structure  to  the  total  damping  of  the  ship  structure  system. 

Betts  and  his  colleagues  conducted  a  survey  of  internal  hull  damping,  and  they 
concluded  that  welding  effects  together  with  stress  concentrations  were  among  the  most 
important  sources  of  hull  damping  in  deformation  modes  [Ref  1].  Figure  1  shows  a 
typical  plot  for  the  damping  properties  of  mild  steel  in  terms  of  the  specific  damping 
strain  energy,  for  a  wide  range  of  stress  amplitudes. 


) 


Amplitude  of  alternating  atresn.OaCMM/m1) 

Figure  1.  Unit  Damping  Energy  vs.  Stress  Amplitude  (from  Ref  7) 


Damping  curves  for  other  metallic  materials  generally  have  similar  forms  to  those 
of  mild  steel.  As  shown  in  Figure  1,  as  stress  amplitude  increases,  damping  energy 
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increases.  The  primary  effect  of  welding  is  from  contraction  of  the  weld  upon  cooling. 
This  contraction  causes  tensile  yield  stresses  in  the  material  parallel  and  adjacent  to  the 
weld.  The  area  of  residual  tensile  yield  stress  at  fillet  welds  usually  extends  for 
approximately  3  to  4  lA  thicknesses  on  either  side  of  the  weld,  leaving  some  10-15%  of 
the  plate  material  at  tensile  yield  stress  [Ref.  1].  Betts  suggests  that  the  residual  stresses 
will  frequently  be  in  the  plastic  region  and  that  vibratory  motion  of  the  weld  will  cause 
plastic  strain,  considerably  increasing  the  damping  of  the  structure.  Imperfect  welds 
where  there  are  gaps  between  the  materials  was  also  believed  to  give  rise  to  dry  friction 
between  the  materials  creating  another  source  of  damping. 

In  thesis  work  done  by  Carey  [Ref.  2]  research  was  conducted  of  the  welding 
effects  on  damping  in  several  beam-stiffened  plates.  At  frequencies  below  500  Hz,  the 
trend  found  was  that  welding  caused  an  increase  in  damping  in  the  plates.  The  results  of 
these  tests  however  were  not  consistent  and  not  precise  enough  to  quantify  the  damping 
attributed  to  weld  effects.  In  a  continuation  of  work  done  by  Carey,  this  thesis  continues 
examination  of  weld  effects  on  damping  in  a  more  complicated  structure  that  contains 
larger  amounts  of  welded  surface,  better  representing  a  ship’s  hull. 

B.  OBJECTIVES 

The  objective  of  this  research  was  to  investigate  the  damping  effects  of  welds  in  a 
beam-stiffened  plate  gaining  insight  on  damping  effects  of  welds  in  a  ship’s  hull 
structure.  A  more  complete  understanding  of  structural  damping  in  ship  structures  may 
improve  computer  simulations  of  ship  vibration  in  ship-shock  trials  and  produce 
opportunities  to  reduce  acoustic  transmission  of  ship  noise. 
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II.  THEORY 


A.  DAMPING  MODEL  ASSUMPTIONS 

As  previously  suggested  in  the  background  section,  the  origins  of  ship  structural 
damping  are  believed  to  lie  in  material  hysteresis  and  dry  friction  at  joints.  Hysteretic 
damping  comes  from  imperfections  in  real  materials  that  cause  energy  to  be  absorbed 
during  cycling,  while  Coulomb  damping  comes  from  the  dry  friction  at  joints. 
Recognizing  the  damping  mechanisms  of  the  system  come  from  a  combination  of 
Hysteretic  and  Coulomb  as  well  as  Viscous  damping,  a  fairly  accurate  mathematical 
representation  can  be  made  by  assuming  only  viscous  damping.  This  assumption  makes 
calculation  of  damping  properties  of  the  model  less  burdensome  while  still  maintaining  a 
good  approximation  to  the  true  total  damping  of  the  model  [Ref.  3]. 

Another  assumption  made  to  make  damping  calculations  less  burdensome  was  to 
assume  the  system  experiences  linear  proportional  damping,  the  damping  of  the  system 
being  linearly  proportional  to  the  mass  and  the  stiffness  of  the  system.  For  a  multi- 
degree-of-freedom  system  such  as  a  ship  hull,  these  assumptions  lead  to  the  following 
damping  equation; 

[C]  =  a[M]  +  /3[K]  (1) 

These  assumptions  fonn  the  basis  for  the  model  by  which  the  damping  values  for 
the  ship  hull  model  were  calculated. 

B.  VISCOUS  DAMPING 

Viscous  damping  can  be  expressed  by  the  following  equation 

Fd  =  cx  (2) 

where  c  is  a  constant  of  proportionality 

The  equation  of  motion  for  a  damped  single  degree  of  freedom  (SDOF)  system 
experiencing  free  vibration  is: 

mx  +  cx  +  kx  =  F(t )  (3) 

Where  m  =  Mass  constant 

c  =  Damping  constant 

k  =  Stiffness  constant 
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x  =  Displacement  of  the  system 
x  =  Velocity  of  the  system 
x  =  Acceleration  of  the  system 
F(t)  =  Force  applied  to  the  system 


The  homogenous  solution  ( F(t)  =  0 )  of  the  differential  equation  of  motion, 

(equation  (3))  is  found  by  assuming  x  =  est .  Substituting  into  the  differential  equation 
yields: 

(ms2  +  cs  +  k)est  =  0  (4) 

which  is  satisfied  for  all  values  of  t  when 

s2+—s  +  —  =  0  (5) 

m  m 

the  two  roots  of  this  characteristic  equation  are: 

s 

Thus  the  general  solution  of  the  differential  equation  is: 

x  =  Ae'1  +  Be2t  (7) 

where  A  and  B  are  constants  evaluated  from  the  initial  conditions  of  displacement  and 
velocity  at  time  zero  (x(O)and  x(0)).  Substituting  equation  (6)  into  equation  (7) 
produces: 


h  c  A2 


1,2 


2m 


2m 


k 

m 


(6) 


x  =  e^c,2m^ 


^  l\J(c/2m)~  —k / m\t  i^j(cl2m)2-k/m\t  ^ 

Ae'  1  +Bex  1 


(8) 


V  ) 

When  the  damping  term  (c/2mf  is  more  than  (A/  in ) ,  the  over  damped  case  will 
occur,  the  exponents  in  (8)  being  real  numbers  and  thus  leading  to  no  oscillations.  If  the 
damping  term  is  less  than  (k!  m ) ,  the  exponents  will  become  imaginary  numbers  leading 


to  the  under  damped  case  and  oscillatory  motion.  When  (cl 2m)  =  kl m  and  the  radical 
is  equal  to  zero  the  case  is  known  as  critical  damping,  cc  and  leads  to  one  oscillation 
cycle 


c 


c 


2mm  =  2y[km 


(9) 
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Any  damping  can  be  expressed  in  terms  of  the  critical  damping  by  a  non- 


dimensional  number  C,  ,  the  damping  ratio,  where  C,  = 


Hence  equation  (6)  becomes 


Sla=(-i±4Fz T)®„  (10) 

and  the  differential  equation  of  motion  (equation  (3))  can  now  be  expressed 

x  +  2£(onx  +  (on  =—  F(t)  (11) 

m 

[Ref  3] 


C.  FORCED  VIBRATION 

To  accurately  calculate  the  damping  in  a  structure,  its  frequency  response  must  be 
determined  in  order  to  obtain  the  natural  frequencies  of  the  structure.  In  any  liner  system, 
there  is  a  direct  liner  relationship  between  the  input  signal  and  the  output  signal;  the  ratio 
of  this  relationship  is  tenned  the  Frequency  Response  Function,  H(co). 

As  previously  stated,  the  equation  of  motion  for  a  single  degree-of-freedom 
system  (1-DOF),  is: 

mx  +  cx  +  kx  =  F(t)  (3) 

To  solve  the  differential  equation,  we  let  the  input  F(t)  =  e‘m' .  The  steady-state  output 
will  be  x  =  H(co)eHOt ,  where  co  is  the  frequency  of  the  applied  force  and  t  is  time  and 
H(co)  is  a  complex  function.  The  steady-state  output  equation  is  differentiated,  to  get 
expressions  for  the  velocity  and  acceleration.  These  expressions  for  displacement, 
velocity  and  acceleration  are  substituted  into  equation  (3)  and  canceling  like  terms  the 
following  is  obtained: 

(-mar  +icco  +  k)H(co)  =  1  (12) 

Then  the  frequency  response  equation  becomes: 

H(co)  =  - - (13) 

k-mco  +ica> 

By  factoring  the  stiffness  from  the  denominator  and  substituting  the  following  equations 
into  equation  (3), 
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m 


(14) 


con  =  Natural  frequency 
£  =  Damping  ratio 
cc  =  Critical  damping  coefficient 
the  frequency  response  can  then  be  written  in  its  classical  form: 


For  a  1-DOF  system,  there  is  one  natural  frequency  and  one  damping  ratio 
associated  with  it.  For  multiple  degrees-of-freedom  (N-DOF),  there  are  as  many  natural 
frequencies  and  damping  ratios  as  there  are  DOF.  Modal  analysis  can  be  used  to  analyze 
a  N-DOF  system.  The  equations  of  motion  for  a  N-DOF  system  in  matrix  form  is: 

[M]{x}  +  [C]{i:}  +  [^]{x}  =  {F}  (16) 

The  mass,  damping,  and  stiffness  matrices  are  n  x  n  matrices,  where  n  is  the  number  of 
DOF  in  the  system.  The  force,  displacement,  velocity  and  acceleration  vectors  are  n  by  1 
in  size.  Each  element  in  the  vectors  corresponds  to  a  DOF  of  the  system.  The  mass  and 
stiffness  matrixes  are  symmetric  and  may  have  some  fonn  of  coupling.  The  first  step  to 
analyze  the  multi-DOF  system  is  to  determine  the  natural  frequencies  and  mode  shapes 
by  analyzing  the  free  response  of  the  system.  The  general  form  of  the  free  response  of 
motion  for  a  N-DOF  system  is  as  follows: 

[M]{x}  +  [if]{x}  =  {0}  (17) 


For  each  of  the  DOF  we  use  x  =  Xem‘  and  equation  (17)  reduces  to: 

Ktj  - M bjO)1  ] { X }  =  0  where  /=  1  ...n  and  /  =  \  ...n  (18) 

The  indices  i  and  j  correspond  to  the  element  locations  in  the  mass  and  stiffness  matrices. 
From  equation  (18),  a  solution  of  the  displacements  is  {X}  =  0  if  the  matrix 

K..-M..co 2  is  invertible.  This  solution,  however,  is  the  trivial  solution.  To  ensure  that 

the  matrix  is  not  invertible,  the  determinant  of  the  matrix  is  forced  to  equal  zero.  By 
forcing  the  matrix  determinant  to  equal  zero,  the  matrix  will  be  singular  and  an  inverse 
matrix  does  not  exist,  thus  a  non-trivial  solution  can  be  found.  The  eigenvalues  found  for 
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the  matrix  are  the  system’s  natural  frequencies  and  the  eigenvectors  found  are  the 
system’s  mode  shapes.  The  mode  shapes  of  the  system  illustrate  how  the  system 
responds  to  an  excitation  at  the  corresponding  natural  frequency.  The  system’s  modal 
matrix  is  formed  by  placing  the  mode  shape  vectors  as  the  columns  of  the  matrix. 

[®]=[WW!-W]  <19> 

®  is  the  modal  matrix  and  (f>  are  the  mode  shape  vectors.  The  modal  matrix  will  have 
the  same  number  of  the  rows  as  there  are  DOF’s.  To  decouple  equation  (16),  we  assume  a 
set  of  modal  coordinates: 

{x}  =  [®]{g},  {x}  =  [®]{g},  {x}  =  [® ]{«)•}  (20) 

next  substitute  equation  (20)  into  equation  (16)  and  multiply  both  sides  by  the  transpose 
of  the  modal  matrix. 

[®rM[®]{»v[*nc][*]{i}+[*n*][®]w=[*r^}  <2i> 

Using  the  orthogonal  properties  of  the  modal  matrix  and  the  symmetric  properties  of  the 
mass,  damping  and  stiffness  matrix  results  in  the  following  equation: 

[mu ]  {<?}  +  [c„  ]  {q}  +  [ku ]  [q]  =  [® f  {F}  (22) 

The  new  modal  mass,  damping,  and  stiffness  matrixes  are  diagonal.  As  a  result,  the 
modal  coordinates  are  decoupled  and  can  be  solved  for  each  Degree  of  Freedom  in  the 
same  manner  as  a  1-DOF  system.  Equation  (22)  can  be  further  simplified  by  multiplying 
both  sides  by  the  inverse  of  the  modal  mass  matrix  to  get: 

=  (23) 

where  Jf|  =[m] 

[Ref.  3  and  4] 

D.  DAMPING  CALCULATION:  HALF-POWER  POINT  METHOD 

The  damping  ratio  can  be  calculated  experimentally  by  using  the  half-power  point 
method.  The  half-power  point  method  detennines  the  damping  ratio  by  examining  the 
sharpness  of  the  resonance  peak.  The  following  equation  is  the  magnitude  of  the 
frequency  response: 
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1- 


+  2^ 


At  resonance,  {&>/&>„  =1)  the  magnitude  of  the  frequency  response  is  Hres  =  —  .  Taking 

2<T 

the  square  of  both  sides  of  equation  (24): 


+  2<r 


+(>-8^)=° 


Solving  for  ( &>/ con )  results  in  the  following  equation: 


(!-2(T2): 


Assuming  that  g  «:  1 ,  the  higher  order  terms  can  be  neglected,  resulting  in  the  following 


equation: 


—  =1±2^ 
co„  , 


Letting  cox  and  co2  correspond  to  the  roots  of  equation  (27)  and  co2  >  o\  ,  equation  (27) 
becomes: 

I 


The  damping  ratio  can  now  be  detennined  by  rearranging  equation  (28): 

_  —  cox  _  /)  ~  f\ 

2  '  2/„ 


Figure  2  illustrates  how  the  half-power  point  method  is  utilized.  [Ref.  1] 
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Amplitude  (dBl 


Frequency  (Hz) 

Figure  2.  Half-power  point  method  [Ref.  1] 


E.  DAMPING  CALCULATION:  COMPLEX  EXPONETIAL  METHOD 


The  Complex  Exponential  Method  (CEM)  is  used  to  extract  modal  parameters  of 
a  system  directly  from  the  system  response  data.  The  Receptance  FRF 
( displacement/ force)  a(co),  of  a  viscously  damped  general  MDOF  system  can  be 
represented  as: 


a 


N  f 

M  =  Z 


-+- 


A 


r=\\JOJ-Sr  JO)-SrJ 


;  where  5,.  =  -©,£•+  My/1- £ 


(30) 


or 


2*  A 

a(co)  =  V - L —  ;  sr  s*r ,  Ar  A*r ,  for  r  >  N 

r= 1  j(0~Sr 

and  Mobility  (velocity/force)  Y{co)  can  be  related  to  a(co) , 


Y(a>)  =  jcoa{co)  (31) 

From  classical  theory,  the  corresponding  Impulse  Response  Function  (IRF),  h(t) 
can  be  obtained  by  taking  the  Inverse  Fourier  Transform  (IFT)  of  the  Receptance  a{co) : 


2  N 

(32) 

r=\ 


The  velocity  form  of  IRF  can  be  expressed  as: 

2  N 

ii{t)  =  YJA,.sres't  (33) 

r= 1 


Then  the  sampled  velocity  data  set  can  be  expressed  as  follows: 
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,h2,....,hq  =  h(0),h(At),h(2At), . ,h(qAt )  (34) 

Using  the  following  simplified  notation, 

es'M  Vr  (35) 

Thus  for  the  jth  sample  data  of  equation  (33)  is  expressed  as: 

2  N  2  N 

h,=hjj& ff  =  =  ZMK  (36) 

r= 1  r= 1 

When  extended  to  the  full  data  set  of  q  samples,  equation  (36)  gives: 

K  —  SlA  S2A  4  h  S2NA2N 

hx  —  VlslAl  +  V2s2A2  4  t  V2Ns2NA2N  (37) 

h2  =  VyxAx  4-  F  s2A2  4  f  V1ns2NA2N 


\  =  V\S\A  +  ^AlA  +  •  •  •  +  V*nS2nAn 

Given  that  the  number  of  sample  points  q  exceeds  4N,  this  equation  can  be  used 
to  set  up  an  eigenvalue  problem  and  the  solution  yields  the  complex  natural  frequencies 
contained  in  the  parameters  V),  V2,  etc.... 

Multiplying  each  equation  in  (37)  by  a  coefficient,  /?.  forms  the  following 

set  of  equations: 

PA  —  PoA  PoA 2  ^  ^  PoAn 

PA  =  pyp  +  PM  +  -  +  Ph2NAN  (38) 

P2h2  =  Ph\  A\  +  Phi  A  + '  • '  +  P2V21nAn 

PA  =  Phq  A  +PAIA+-+ Phi  A, 

Adding  all  equation  in  (38)  results  in: 

q  2  N  q 

£M=Z(4X/A)  (39) 

1=0  7=1  i= 0 

The  coefficients  /?.  are  taken  to  be  the  coefficients  in  the  equation, 

p0  +  py  +  py1  +  py  A-. +pqv“=  o  (40) 

for  which  the  roots  are  V1 ,  V2 , . . . ,  V  . 
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The  values  of  the  ft  coefficients  are  sought  to  detennine  the  roots  of  (40),  the 
values  of  Vr  and  thus  the  system  natural  frequencies.  With  q  being  the  number  of 

degrees  of  freedom  of  the  system  model,  we  can  set  these  two  parameters  to  the  same 
value,  i.e.  let  q=2N. 

The  equation  (40)  can  now  be  expressed  as: 


2N 


'YjPjK-  =  0  ;  for  1=1,  2N 

7=0 

and  thus  every  term  on  the  right-hand  side  of  equation  (39)  is  zero. 

2N 

£M=» 

1=0 

rearranging  equation  (42),  by  setting  J32N  =  1 

in-\ 

X  PA  =  ~Kn 

1=0 


(41) 


(42) 


(43) 


o 

1 _ 

K 

••  h 

n2N-\ 

'  00  ' 

f 

v _ 

K 

•-sT  ... 

K  • 

...  .i5- 

% 

< 

Pi 

>  =  —  < 

h 

n2N+\ 

h 

[_n2N-\ 

h 

n2N 

h 

n2N+\ 

••  h  1 

r,4N- 2  J 

02N-1 , 

h 

ln4N-\  J 

Repeating  the  process  from  (34)  to  (43)  using  different  sets  of  IRF  data  points  and 
further  choosing  new  data  sets  that  overlaps,  successive  applications  of  this  procedure 
lead  to  a  full  set  of  2N  equations: 


(44) 


or 

' h 1  {/?}  =-(£  1  (45) 

-  J2Nx2N  t^t2JVxl  1  )2Nxl 

The  unknown  coefficients  {/?}  can  be  found  from  equation  (45).  The  values  of 
V1,V2,...,V2N  can  now  be  detennined  using  equation  (40)  and  subsequently  the  system 
natural  frequencies  can  be  found  using  the  relationship. 

Vr=esM  (35) 

Using  equation  (37),  corresponding  modal  constants,  A1,A2,...,A2N  can  be 
calculated.  This  may  be  written  as: 
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1 

1 

1 

1 

o 

J 

K 

V2 

V2 

V 

r  2  N 

A-2&2 

K 

v2 

V2 

v; 

y2 

v  2N 

5 

A3S3 

>  =  < 

K 

V™-1 

yiN- 1 

v  2 

y2N-\ 

...  v2N~l  1 

v  2  N  J 

^A2NS2N> 

h 

L  2JV— 1  J 

or 

[V]{A}  =  {h)  (47) 

[Ref.  5] 


F.  VERIFICATION  OF  EXTRACTED  MODAL  PARAMETERS 


The  modal  parameters  calculated  from  the  above  procedure  can  be  verified  by 
comparing  synthesized  time  histories  to  the  originally  measured  time  histories.  From 
equations  (30)  and  (31), 


jo)Ar  jcoAi 


+  - 


"W® 


-s„ 


jo)-sr) 


Sr=-OrCr+jG)r^-Cr 


(48) 


or 

Y(a>)  =  V  ^  s*r ,  Ar  ->  4  ,  for  r  >  N  (49) 

r=l  j(0~Sr 

Mobility  can  be  calculated  and  by  inverse  FFT,  taking  the  real  part  of  the  results, 
synthesized  IRF  can  be  calculated  and  compared  to  original  time  histories.  [Ref.  5] 

j/zj  =  Real  of  IFFT  (A/  Y)  (50) 

G.  CALCULATION  OF  RAYLEIGH  DAMPING 

Using  Rayleigh  damping  representation,  the  damping  matrix  can  be  represented 
as: 

[C]  =  a[M]  +  p[K]  (1) 

by  using  the  mass  normalized  modal  matrix  [^] ,  the  damping  matrix  may  we  rewritten: 

Wr[c]M=[2».cU (5i) 

By  using  equation  (5 1)  for  all  2N  modes,  following  2N  equations  can  be  set. 
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a  +  fio)]  =  2  a£x 

a  +  pm;=  l(o2C2  ^ 

oc  +  P<x>2N  —  2o)2NC2N 


or 


(53) 


From  equation  (53),  the  two  parameters  a  and  ft  are  calculated  in  a  least  squares 


solution  by  matrix  pseudo-inversion.  [Ref.  5] 


(MLRLjVl 


T 

2x2  N 


{Z)lNxl 


(54) 
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III.  EXPERIMENTAL  SETUP 


A.  EQUIPMENT  DESCRIPTION 

The  equipment  used  during  testing  was  as  follows: 

•  8.5  ft  x  4  ft  “Hull  model”  panel 

•  8.5  ft  x  1.5  ft  flat  panel 

•  Hewlett  Packard  3  5  62 A  Dynamic  Signal  Analyzer 

•  PCB  Piezotronics  Modally  Tuned  Impact  Hammer 

•  PCB  Piezotronics  model  336C04  Accelerometers 

•  PCB  Piezotronics  model  208C01  Force  Sensor 

•  PCB  Piezotronics  model  483B07  Power  Unit  Amplifier 

•  MB  Dynamics  model  A  SS250VCF  Amplifier 

•  MB  Dynamics  Modal  50  Exciter 

•  Modified  PC  computer  with  National  Instruments  “Lab VIEW”  software 

Experimental  testing  was  conducted  on  a  stiffened  metal  panel  matching  a  Vi  scale 
section  of  generic  warship  hull  structure.  The  model  plate  was  4  ft  by  8.5  ft,  1 1  gauge  A- 
36  steel  and  was  stiffened  length-wise  and  around  the  plate  edges  by  1  inch  high  17 
gauge  A-36  steel  strips  with  3  inch  spacing  from  the  edge  of  the  plate  and  6  inch  periodic 
spacing  between  stiffeners.  Transverse  stiffeners  were  1.75  inch  high  14  gauge  A-36 
steel  strips  with  3  inch  spacing  from  the  edge  and  24  inch  periodic  spacing.  Stiffeners 
were  double  filet  welded  to  the  plate.  The  hull  model  was  freely  suspended  from  an  H- 
frame  by  a  series  of  nylon  straps  hooked  to  steel  rings  welded  to  the  structure. 
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3” 


Figure  4.  Photo  of  panel 

A  section  of  flat  panel  1.5  ft  by  8.5  ft,  of  4  gauge  A-36  steel  was  tested  to  be  used 
as  a  reference  to  compare  with  the  behavior  of  the  multi-weld  “hull  panel”.  The  flat  plate 
was  similarly  suspended  by  a  series  of  bungee  cords  to  simulate  a  free  system 
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Figure  5.  Schematic  of  Flat  Panel  model  and  photo  of  panel 

The  Hewlett  Packard  3562A  Dynamic  Signal  Analyzer  (DSA)  was  used  to 
compute  the  frequency  response  of  the  panels  and  as  a  signal  generator.  The  DSA  is  a 
dual-channel,  fast  Fourier  transform-based  network,  spectrum  and  waveform  analyzer 
that  provides  analysis  capabilities  in  both  the  time  and  frequency  domains.  Two  input 
channel  and  a  built  -in  signal  source  (noise  and  sine  signals)  can  be  used  to  perform 
spectrum  analysis.  The  DSA  has  a  frequency  resolution  of  25.6  pHz  allowing  the  user  to 
obtain  highly  accurate,  high-resolution  plots  of  the  frequency  responses  of  the  mechanical 
system.  Single  channel  accuracy  is  ±0.15  dB  with  80  dB  of  dynamic  range  [Ref.  6]. 

The  Piezotronics  Modally  Tuned  Impact  Hammer  was  the  series  086B03  model. 
The  modally  tuned  hammer  is  capable  of  exciting  structural  resonances  up  to  10  KHz  and 
has  a  sensitivity  of  10  mV/lbf.  The  hammer  converts  a  transfer  force  into  an  electrical 
signal  which  was  then  used  by  the  DSA  as  the  input  response.  The  input  signal  was 
obtained  from  a  force  transducer  on  the  head  of  the  impact  hammer. 

The  Piezotronics  model  336C04  Accelerometers  are  hermetically  sealed,  shear 
structured  ICP  accelerometers.  The  accelerometers  provide  a  10  mV/g  output  over  a 
frequency  range  from  1  to20,000  Hz  (±10%). 

PCB  Piezotronics  model  208C01  Force  Sensor  was  used  to  measure  exciter  force 
input  into  the  test  panel.  The  force  sensor  is  an  ICP  Quartz  Force  Sensor  and  measures  in 
a  range  from  10  lbf  tension  to  10  lbf  compression  with  a  sensitivity  of  502.5  mV/lbf  [Ref. 
7  and  Ref  8]. 

The  PCB  Piezotronics  model  483B07  Power  Unit  Amplifier  is  a  12  channel  amplifier 
that  used  to  amplify  the  signals  from  the  force  sensor  and  Accelerometers. 

The  MB  Dynamics  model  A  SS250VCF  Amplifier  was  the  used  to  power  the  Modal 
50  Exciter  and  amplify  the  input  signal  into  the  Exciter. 
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The  MB  Dynamics  Modal  50  Exciter  is  a  lightweight  pennanent  magnet  shaker  that 
provides  50  pounds  dynamic  force  output.  The  force  output  from  the  Exciter  was 
transmitted  to  the  panel  through  a  steel  rod  threaded  into  the  PCB  Piezotronics  force 
sensor  that  was  bolted  to  the  panel. 

A  Personal  Computer  with  a  multi-channel  input  signal  circuit  board  was  used  to  take 
the  signals  from  the  force  sensor  and  accelerometers  via  the  PCB  Piezotronics  Power 
Unit  Amplifier.  National  Instruments  “LabVIEW”  software  was  used  to  pair  the  force 
sensor  with  each  of  the  accelerometers  and  making  fast  Fourier  transforms  on  the  pairs  to 
measure  the  magnitude/phase  and  real/imaginary  frequency  response  and  coherence  of 
the  single  input  multi-output  system. 

B.  EQUIPMENT  SETUP 
1.  Impact  Testing 

For  impact  testing,  the  PCB  modally  tuned  impact  hammer  was  used  to  excite  the 
model  while  accelerometers  (PCB  Piezotronics  model  336C04)  measured  the  structure’s 
response.  The  impact  hammer  test  was  used  to  determine  the  frequency  response  of  the 
model.  The  model  was  impacted  at  the  excitation  point  and  accelerometer  readings  were 
taken  at  a  variety  of  positions  on  the  model.  The  input  and  output  signals  were  amplified 
by  the  PCB  Power  unit  with  input  gain  set  at  1  and  output  gain  set  at  10.  The  input  and 
output  signals  were  analyzed  using  the  Dynamic  Signal  Analyzer  (DSA),  where  the 
frequency  response  and  coherency  were  measured.  The  force-exponential  window  was 
used  during  the  impact  hammer  test  and  10  stable  mean  averages  were  used  to  obtain  the 
frequency  response.  The  time  record  was  set  to  trigger  when  the  input  signal  from  the 
hammer  reached  0.5  volts.  Figure  6  shows  the  equipment  setup  for  impact  hammer 
testing. 
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Output  signal 


Figure  6.  Equipment  setup  for  Impact  Hammer  Testing 


2.  Exciter  Testing 

The  Dynamic  Signal  Analyzer  provided  a  1  Volt  random  noise  input  signal 
throughout  the  range  of  0-250  Hz  via  the  MB  Dynamics  amplifier  to  the  MB  Dynamics 
Modal  50  Exciter  to  excite  the  panel.  The  input  and  output  signals  were  amplified  by  the 
PCB  Piezotronics  Power  Unit  Amplifier  with  gain  set  to  1  for  all  channels.  The  input 
signal  was  measured  by  the  PCB  Piezotronics  force  sensor  bolted  to  the  panel  to  which 
the  exciter’s  steel  rod  was  threaded  into.  PCB  Piezotronics  accelerometers  were  used  to 
measure  the  model  response  throughout  the  structure.  The  signals  of  the  input  force 
sensor  and  the  one  output  accelerometer  were  fed  to  the  Digital  Signal  Analyzer  for 
frequency  response,  and  coherency  measurement  and  the  input  signal  and  all  5  output 
signals  were  fed  to  the  PC  where  the  National  Instruments  “Lab View”  software  to 
measure  the  magnitude/phase  and  real/imaginary  frequency  response  and  coherence.  10 
stable  mean  averages  were  used  to  obtain  the  model’s  frequency  response  and  frequency 
resolution  was  set  to  0.05  Hz.  Results  were  saved  as  text  files  for  further  post-processing. 
Figure  7  shows  a  schematic  of  the  experimental  setup. 
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Accelerometers 

Figure  7.  Equipment  setup  for  Exciter  Testing 


C.  FINITE  ELEMENT  MODEL  SIMULATION 

In  an  attempt  to  relate  damping  with  the  mode  shape  and  natural  frequency  a 
finite  element  model  was  created  for  hull  panel  and  flat  plate.  The  panel  was  modeled 
using  22092  quad  4  shell  elements  utilizing  the  MSC  Patran/Nastran  computer  modeling 
system.  No  attempt  was  made  to  simulate  the  actual  weld  of  the  structure,  stiffener  and 
panel  surfaces  were  simply  merged  together  using  the  equivalence  command  in  Patran. 
Similarly,  the  flat  plate  was  modeled  using  29376  quad  4  shell  elements. 
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IV.  EXPERIMENTAL  SETUP  AND  RESULTS 


A.  IMPACT  HAMMER  TESTING 

A  variety  of  tips  of  varying  hardness  were  used  with  the  impact  hammer  to  find 
which  tip  produced  the  best  coherence.  During  testing,  an  acceptable  level  of  coherency 
was  not  able  to  be  achieved  using  any  of  the  various  tips  for  the  impact  hammer  and 
impact  hammer  testing  was  abandoned  in  favor  of  shaker  testing  which  produced  much 
better  levels  of  coherence. 

B.  EXCITER  TESTING 

1.  Half-Power  Point  Method 

The  Half-Power  point  method  was  used  to  make  preliminary  tests  to  find  the 
effect  damping  across  the  welded  panel.  For  this  initial  testing  the  force  gage  and  a 
single  accelerometer  were  used  as  the  input  and  output  respectively.  Figure  7  displays 
the  placement  of  the  force  gage  and  the  various  placements  of  the  accelerometer  across 
the  panel  for  this  preliminary  testing. 


Figure  8.  Exciter  (*)  and  Accelerometer  (letters)  Placement  for  Half- 

Power  Point  Testing 

A  1  Volt  random  noise  signal  from  the  DSA  was  used  as  the  input  signal.  The 
testing  frequency  range  was  from  0  to  250  Hz  and  was  averaged  15  times.  Figure  8  is  a 
sample  frequency  plot  of  the  hull  panel  response. 
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Through  testing  a  resonant  peak  was  found  as  approximately  113  Hz  throughout 
the  panel.  Using  a  zoom  measurement  of  the  113  Hz  frequency  range,  a  more  precise 
measurement  of  the  resonance  frequency  fn ,  and  the  frequencies  f\  and  f2 ,  3  decibels 

below  resonance  before  and  after  resonance  respectively  were  found.  As  discussed  in  the 
theory  section  previously,  the  damping  ratio  can  now  be  calculated  using  equation  (29) 


„  _  (O2  _  ft  f\ 

“  2  fn 

Table  1.  Half-Power  Point  Testing  Results 


f2  =113.137  A 

f2  =  113.05  H 

/2  =112.962 
f  =112.817 

fn  =113-11 

O 

/2  =113.1 

f  =112.97  0  00146 

fn  =113-3 

/  =112.9  f  ~  0.00133 

fn  =113-2 

=0.00130 

f=  112.93  ^=00°153 

/„  =113.275 

f2  =113.162  B 

f2  =  113.137  1 

/2  =113.112 

P 

/2  =113.15  W 

/  =113  f  —  0.00137 

y;  =  112.987  1^=0.00130 

/j  =112.968 

=0.00133 

/j  =113  =0.0013 

fn  =113-31 

f  =113.281 

«/  71 

f  =113.268 

=113.312 

f2  =113.212  C 

f2  =  113.137  J 

/2  =113.075 

Q 

/2  =113.212  X 

fl=  113.03  f  —  0.00158 

/  =  112.975  Q  —  0.00144 

f  =112-9 

=0.00155 

/  =  113.037  tT  =  0.00199 

fn  =113.387 

fn  =  II3-3 

/„  =113-25 

/n  =113.487 

f2  =113.125  ° 

f2  =  113.162  K 

/2  =113.125 

R 

/2  =113.175  Y 

fx  =112.95  ^=00°133 

/  =113.06  =  0.00117 

y  =112.969 

^  =0.00138 

/j  =  113.018  <y=°-°°i38 

fn  =  113-25 

f  =113.325 

./  « 

f  =113.282 

f  =113.331 

•/  « 

f2=  113.175  E 

/2  =  113.062  L 

/2  =112.9 

S 

/2  =113.175  Z 

fx  =113.012  f  —  0.00133 

/  =  112.906  =0.00134 

f  =112.73 

=0.00147 

/  =  113.025  4"  =  0.00133 

fn  =113-312 

fn  =113.21 

f  =113.062 

J  n 

f  =113.325 

J  n 

/2  =113.162  F 

/2  =113.175  M 

f2  =113.15 

T 

f2  =113-137  AA 

f-  113.06  =0.00119 

/  =113.018  =  0.00136 

f  =113.018 

^  =0.00125 

f=  112.975  ^=0.00144 

/„  =11333 

fn  =113.325 

fn  =113.3 

/„  =113-3 

/2  =113.15  G 

fl  =113-1  N 

/2  =113.062 

U 

/2  =113.137  BB 

/j  =  112.98  f  —  0.00146 

/  =112.94  =  0.00137 

f  =112-9 

=0.00137 

/  =  112.98  4"  =  0.00141 

fn  =113-31 

fn  =113-25 

/„  =113-21 

/„  =113-3 

Table  1  shows  the  results  of  half-power  point  testing.  Of  particular  note  is  the 
over-all  low  level  of  damping  throughout  the  structure,  all  positions  displaying  less  the 
0.2%  damping. 
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2.  Modal  Parameter  Extraction 

For  more  robust  analysis  of  the  frequency  dependence  of  damping  as  well  as 
position  in  the  panel  (and  thus  the  wave  propagation  through  welds)  the  modal  parameter 
extraction  technique  was  used.  Similarly  to  the  half-power  point  method,  a  1  Volt 
random  noise  signal  from  the  DSA  was  used  as  the  input  signal  to  the  exciter.  Two 
frequency  ranges  were  tested,  0  to  250  Hz  for  a  comprehensive  look  at  vibration  through 
the  frequency  range  of  interest  and  0  to  20  Hz  for  a  more  studied  examination  of  the 
lower  frequencies.  Frequency  resolution  of  0.05  Hz  and  10  stable  mean  averages  were 
used  for  both  the  broad  band  and  narrow  band  tests.  Using  the  National  Instruments 
“Lab View”  software  on  the  PC,  five  accelerometers  at  a  time  (output  signals)  were 
processed  with  the  input  signal  (force  sensor)  to  find  fast  Fourier  transforms,  calculating 
magnitude/phase  and  real/imaginary  frequency  response  and  coherence  for  each 
input/output  pair.  The  processed  fast  Fourier  transforms  results  were  saved  as  text  files 
for  post-processing  and  analysis  of  the  FFT  data. 

The  modal  parameter  extraction  and  verification  of  extracted  modal  parameters 
techniques  articulated  in  the  theory  section  were  used  in  a  Fortan  code  to  calculate  the 
modal  parameters  of  the  position  and  regenerated  the  FFT  curve  using  Microsoft 
Developer  Studio’s  Fortran  PowerStation  4.0.  Mircocal  Origin  6.0  was  used  to 
comparing  the  FFT  curve  and  regenerated  FFT  curves.  The  systems  most  resonant 
modes  and  their  associated  damping  were  found  and  then  saved  to  an  Excel  spreadsheet 
for  positional  comparison  of  frequency  dependant  damping.  Figures  9  and  10  show  an 
examples  of  resonant  modes  that  were  chosen  for  positional  comparison. 
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Magnitude 


Posit  37  -  20  Hz 


Frequency  (Hz) 


Figure  9.  Sample  Frequency  response  and  Regenerated  Frequency  response. 
Circled  are  strongly  resonant  modes  that  would  be  chosen  for  positional 

comparison. 
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Posit  37  -  250  Hz 


Figure  10.  Sample  Frequency  response  and  Regenerated  Frequency  response. 
Circled  are  strongly  resonant  modes  that  would  be  chosen  for  positional 
comparison. 

Rayleigh  damping  a  and  (3  for  each  position  and  the  resulting  curve  fit  was 
calculated  as  articulated  in  the  Theory  section  by  using  the  modal  frequencies  and 
damping  ratios  found.  MATLAB  code  was  used  to  use  calculate  these  damping 
coefficients  and  to  plot  the  curve  fit  compared  to  measured  damping. 

a.  Horizontal  Position  Examination 

Positional  relation  of  damping  compared  horizontally  from  the  exciter  is 
shown  in  Figure  1 1 
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Figure  11.  Horizontal  Positions  Compared 


The  testing  positions  along  the  longitudinal  stiffener  (Positions  7,  22,  37,  and 
52)  and  those  off  stiffener  (Positions  6,  21,  36,  and  51)  were  compared  separately  as 
show  on  figures  12  and  13  and  then  combined  in  figure  14. 


Horizontal  Weld  Damping 


Frequency  (Hz) 

♦  Posit  7  ■  Posit  22  Posit  37  Posit  52 


Figure  12.  Horizontal  Weld  Damping 
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Horizontal  Off-Weld  Damping 


Frequency  (Hz) 

♦  Posit  6  r  Posit  21  Posit  36  Posit  51 


Figure  13.  Horizontal  Off-Weld  Damping 


Frequency  (Hz) 

♦  Posit  6  ■  Posit  21  Posit  36  Posit  51  x  Posit  7  •  Posit  22 

+  Posit  37  -  Posit  52 

Figure  14.  Total  Horizontal  Damping 

As  seen  from  Figures  12,  13  and  14  there  is  little  difference  in  damping  on  or 
off  the  weld  and  little  difference  in  damping  in  relation  to  testing  position.  Damping  as  a 
general  trend  however,  is  frequency  dependant  displaying  an  exponential  decay  with 
frequency. 

The  calculated  Rayleigh  damping  a  and  f5  for  each  position  are  shown  in 
Table  2  and  the  resulting  curve  fit  compared  to  the  original  data  points  are  shown  in 
Figures  15,  16,  17,  and  18. 
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Table  2.  Rayleigh  Damping  Results  for  Horizontal  Damping 


Position 

Alpha 

Beta 

6 

0.5818 

2.6823E-6 

7 

0.5921 

2.6804E-6 

21 

0.7358 

2.7782E-6 

22 

0.7163 

2.6219E-6 

36 

0.5165 

2.4298E-6 

37 

0.6539 

3.072  IE-6 

51 

0.6724 

2.4508E-6 

52 

0.6947 

2.8037E-6 

The  mean  a  calculated  for  horizontal  positions  was  0.6454,  with  a  standard 
deviation  of  0.0755  and  the  mean (3  calculated  for  horizontal  positions  was  2.6899E-6 
with  a  standard  deviation  of  2.0577E-7. 


Posit  6 


Figure  15.  Positional  Damping  Rayleigh  Damping  Curve-Fit 
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Posit  51 


Figure  18.  Positional  Damping  Rayleigh  Damping  Curve-Fit 


b.  Vertical  Position  Examination 

Positional  relation  of  damping  compared  vertically  near  the  exciter  is  shown 
in  figure  19. 


Figure  19.  Vertical  Positions  Near  Exciter  Compared 
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The  vertical  testing  positions  near  the  exciter  (Positions  21,  22,  23,  24,  25,  and 
26)  were  compared  on  figure  20. 


Vertical  Damping  Near  Exciter 


Frequency  (Hz) 


■  Posit  21  •  Posit  22  ♦  Posit  23  Posit  24  Posit  25  x  Posit  26 

Figure  20.  Vertical  Damping  Near  Exciter 

Similarly  to  the  Horizontal  position  examination,  figure  20  shows  there  is 
little  difference  in  relation  to  testing  position  and  damping  generally  displayed  an 
exponential  decay  with  frequency. 

The  calculated  Rayleigh  damping  a  and  /?  for  each  position  are  shown  in 
table  3  and  the  resulting  curve  fit  compared  to  the  original  data  points  are  shown  in 
Figures  21,  22,  and  23. 


Table  3.  Rayleigh  Damping  Results  for  Vertical  Damping  Near  Exciter 


Position 

Alpha 

Beta 

21 

0.73580 

2.7782E-6 

22 

0.7163 

2.6219E-6 

23 

0.7376 

2.4544E-6 

24 

0.7332 

2.7233E-6 

25 

0.7361 

2.3462E-6 

26 

0.5428 

2.4197E-6 
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The  mean  a  calculated  for  vertical  positions  near  the  exciter  was  0.7003  , 
with  a  standard  deviation  of  0.0776  and  the  mean /?  calculated  for  vertical  positions  near 
the  exciter  was  2.5573E-6  with  a  standard  deviation  of  1.7586E-7. 
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Positional  relation  of  damping  compared  vertically  far  from  the  exciter  is 
shown  in  figure  24 


Figure  24.  Vertical  Positions  Far  from  the  Exciter  Compared 

The  vertical  testing  positions  far  from  the  exciter  (Positions  51,  52,  53,  54,  55, 

and  56)  are  compared  on  figure  25. 


Vertical  Damping  far  from  Exciter 


Frequency  (Hz) 


Posit  51  ■  Posit  52  ♦  Posit  53  ■  Posit  54  Posit  55  x  Posit  56 


Figure  25.  Vertical  Damping  far  from  Exciter 
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Again,  similar  to  the  horizontal  position  examination  and  the  vertical  positions 
near  the  exciter,  figure  25  shows  there  is  little  difference  in  relation  to  testing  position 
and  damping  generally  displayed  an  exponential  decay  with  frequency. 

The  calculated  Rayleigh  damping  a  and  /?  for  each  position  are  shown  in 
table  4  and  the  resulting  curve  fit  compared  to  the  original  data  points  are  shown  in 
Figures  26,  27,  and  28. 


Table  4.  Rayleigh  Damping  Results  for  Vertical 


Damping  Far  From  1 

Exciter 

Position 

Alpha 

Beta 

51 

0.6724 

2.4508E-6 

52 

0.6947 

2.8037E-6 

53 

0.7660 

2.4185E-6 

54 

0.7925 

2.2097E-6 

55 

0.8263 

2.1606E-6 

56 

0.5483 

2.8479E-6 

The  mean  a  calculated  for  vertical  positions  near  the  exciter  was  0.7167, 
with  a  standard  deviation  of  0. 1010  and  the  mean /?  calculated  for  vertical  positions  near 
the  exciter  was  2.4819E-6  with  a  standard  deviation  of  2.8978E-7. 
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Posit  51 


0  50  100  150  200  250 

Frequency  (Hz) 


Figure  26.  Positional  Damping  Rayleigh  Damping  Curve-Fit 


Posit  53 


0_ i  i  1  ' I  i  ■*-  *  -T- T' 

0  50  100  150  200  250 

Frequency  (Hz) 


Figure  27.  Positional  Damping  Rayleigh  Damping  Curve-Fit 
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Posit  55 

0.03  , - , - , - 


Frequency  (Hz) 
Posit  56 


Figure  28.  Positional  Damping  Rayleigh  Damping  Curve-Fit 

The  overall  Rayleigh  damping  throughout  the  structure  comparing  all 
positions  found  a  mean  a  value  of  0.6779,  with  a  standard  deviation  of  0.0953  and  a 
mean / 3  value  of  2.5687E-6  with  a  standard  deviation  of  2.4747E-7. 


C.  FINITE  ELEMENT  MODELING 

Using  the  MSC  Patran/Nastran  finite  element  modeling  program,  a  finite  element 
model  was  developed  for  the  panel.  The  attempt  was  to  observe  trends  that  would  relate 
the  mode  shape  of  the  panel  with  damping  and  natural  frequency.  Table  5  shows  the 
natural  frequencies  and  mode  shape  numbers  obtained  from  the  finite  element  model. 


Table  5.  FEM  Calculated  Natural  Frequencies 


Mode 

Frequency  (Hz) 

Mode 

Frequency  (Hz) 

1 

5.304 

18 

146.34 

2 

11.988 

19 

148.58 

3 

21.096 

20 

156.36 

4 

32.2 

21 

159.06 

5 

47.129 

22 

173.78 

6 

60.35 

23 

178.34 

7 

66.92 

24 

181.77 

37 


8 

74.14 

25 

185.28 

9 

82.477 

26 

186.47 

10 

91.227 

27 

198.54 

11 

93.094 

28 

204.16 

12 

111.56 

29 

209.2 

13 

122.25 

30 

211.3 

14 

124.24 

31 

223.2 

15 

129.95 

32 

225.1 

16 

135.15 

33 

230.54 

17 

143.79 

34 

241.91 

Table  6  shows  FEM  natural  frequencies  and  modes  that  matched  highly  resonant 
experimental  natural  frequencies. 


Table  6.  Comparison  of  FEM  and  Experimental  Natural  Frequencies 


FEM  Mode 
Shape  # 

FEM  Natural 
Frequency  (Hz) 

Experimental  Natural 
Frequency  (Hz) 

2 

11.988 

11.3-11.4 

3 

21.096 

24.4-24.6 

4 

32.2 

36.3 

7 

66.92 

66.3-66.6 

12 

111.56 

112.9-113.3 

24 

181.77 

183.0 

28 

204.16 

205.5-206.1 

29 

209.2 

208.5-208.9 

33 

230.54 

232.3-232.9 

Figures  29  to  33  show  examples  of  mode  shapes  found  by  the  FEM  model  that 
matched  highly  resonant  experimental  modes.  Flexural  (transverse  and  longitudinal)  and 
torsional  mode  shapes  as  well  as  combination  mode  shapes  were  all  found  to  match  with 
highly  resonant  experimental  natural  frequencies.  No  type  of  mode  shape  was  found  to 
be  preferential  in  highly  resonant  experimental  natural  frequencies. 
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Figure  29.  FEM  mode  shape  with  natural  frequency  of  11.988  Hz 


Figure  30.  FEM  mode  shape  with  natural  frequency  of  21.096  Hz 


Figure  31.  FEM  mode  shape  with  natural  frequency  of  32.2  Hz 
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Figure  32.  FEM  mode  shape  with  natural  frequency  of  66.92  Hz 


Figure  33.  FEM  mode  shape  with  natural  frequency  of  209.2  Hz 


D.  FLAT  PLATE  COMPARISON 

For  reference  to  compare  against  the  welded  panel,  the  flat  panel  previously 
described  in  the  experimental  setup  section  was  subjected  to  the  same  exciter  testing  and 
modal  parameter  extraction  as  the  plate  stiffened  panel.  Figure  29  displays  the  placement 
of  the  force  gage  and  accelerometers  across  the  panel  for  this  testing. 
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*  1  2  3  4  5 


Figure  34.  Exciter  and  Accelerometer  Placement  for  Flat  Panel  Testing 


The  testing  positions  along  the  flat  plate  are  compared  in  figure  30. 
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Figure  35.  Flat  Panel  Damping 

As  with  the  plate  stiffened  panel,  damping  there  appears  to  be  little  difference  in 
relation  to  testing  position  on  the  flat  plate  and  damping  generally  displayed  an 
exponential  decay  with  frequency.  The  difference  between  the  flat  plate  and  plate 
stiffened  panel  is  the  overall  magnitude  of  damping  for  the  flat  plate  is  lower  than  the 
plate  stiffened  panel  especially  in  lower  frequencies  (less  the  50  Hz). 
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The  calculated  Rayleigh  damping  a  and  [3  for  each  position  are  shown  in  table  7 
and  the  resulting  curve  fit  compared  to  the  original  data  points  are  shown  in  Figures  30, 
31,  and  32. 


Table  7.  Rayleigh  Damping  Results  for  Flat  Plate 


Position 

Alpha 

Beta 

1 

0.3797 

1.0395E-6 

2 

0.3694 

1.0657E-6 

3 

0.3720 

1.0353E-6 

4 

0.3734 

1.0320E-6 

5 

0.3628 

1.0215E-6 

The  mean  a  and  (3  calculated  for  horizontal  positions  were  0.3715  and  1.0388E- 
6  respectively  and  standard  deviations  of  0.0062  and  1.645  IE-8.  This  agrees  with  the 
damping  curves  of  figure  30,  a  being  about  55%  the  value  of  a  compared  to  the 
stiffened  plate  and  f3  being  about  40%  the  value  of  [3  compared  to  the  stiffened  pate. 


Figure  36.  Positional  Damping  Rayleigh  Damping  Curve-Fit 
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Flat  Panel  Posit  3 


Frequency  (Hz) 
Flat  Panel  Posit  4 


Frequency  (Hz) 


Figure  37.  Positional  Damping  Rayleigh  Damping  Curve-Fit 


Flat  Panel  Posit  5 


Figure  38.  Positional  Damping  Rayleigh  Damping  Curve-Fit 
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As  with  the  plate  stiffened  panel,  a  finite  element  model  was  developed  for  the 
flat  plate  in  an  attempt  to  observe  trends  relating  the  mode  shape  of  the  panel  with 
damping  and  natural  frequency.  Table  8  shows  the  natural  frequencies  and  mode  shape 
numbers  obtained  from  the  finite  element  model. 


Table  8.  FEM  Calculated  Flat  Plate  Natural  Frequencies 


Mode 

Frequency  (Hz) 

Mode 

Frequency  (Hz) 

1 

4.913 

13 

124.25 

2 

13.596 

14 

152.84 

3 

16.821 

15 

158.16 

4 

26.772 

16 

164.4 

5 

34.338 

17 

169.77 

6 

44.436 

18 

178.29 

7 

53.20 

19 

185.94 

8 

66.603 

20 

196.43 

9 

73.988 

21 

201.41 

10 

93.247 

22 

219.41 

11 

97.22 

23 

222.93 

12 

123.38 

24 

243.22 

Table  9  shows  FEM  natural  frequencies  and  modes  that  matched  highly  resonant 
experimental  modes. 


Table  9.  Comparison  of  Flat  Plate  FEM  and  Experimental  Natural  Frequencies 


FEM  Mode 
Shape  # 

FEM  Natural 
Frequency  (Hz) 

Experimental  Natural 
Frequency  (Hz) 

1 

4.913 

4.92 

2 

13.596 

13.45 

4 

26.772 

26.55 

6 

44.436 

44.07 

8 

66.603 

66.13 

10 

93.247 

92.71 

15 

158.16 

157.48 

16 

164.4 

164.49 

17 

169.77 

169.67 

18 

178.29 

177.83 

20 

196.43 

195.62 

21 

201.41 

200.55 

22 

219.41 

219.27 

24 

243.22 

242.26 
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The  modes  and  frequencies  of  the  FEM  model  closely  matched  all  resonant 
modes  of  the  flat  plate  found  experimentally  and  only  a  few  “false”  modes  were  created 
by  the  FEM  model.  This  result  confirmed  the  veracity  of  the  FEM  model  but  added  no 
further  insights  into  the  relationship  between  welds  and  damping. 


45 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


46 


V.  CONCLUSIONS  AND  RECOMEND ATIONS 


The  results  of  this  damping  study  indicate  that  welds  do  increase  the  damping  of 
structure  and  that  Rayleigh  damping  is  and  appropriate  model  for  its  analysis.  The  study 
results  indicate  that  damping  is  relatively  high  in  the  low  frequency  domain  while 
decreasing  exponentially  in  the  higher  frequencies.  The  damping  from  the  welds 
increased  the  damping  of  the  overall  system  but  did  not  show  any  positional  damping. 
Damping  values  found  near  the  point  of  excitation  were  similar  to  damping  values  of 
positions  separated  from  the  point  of  excitation  by  multiple  welds. 

This  level  of  damping  increase  found  compare  to  a  bare,  flat  plate  however  is 
extremely  small,  with  average  Rayleigh  damping  coefficients  of  a  =  0.6779  and 
/?  =  2.5687E-6  found.  From  DDG-53  ship  shock  trial  data,  Rayleigh  damping 
coefficients  of  a  =  19.2  and  /?  =  2.09E-6  were  estimated  for  ship  structural  damping. 
[Ref.  9]  The  minor  increase  in  damping  found  due  to  welds  in  this  study  is  unlikely  to  be 
a  significant  contributor  to  overall  this  ship  structure  damping  found. 

Recommendations  for  future  work  into  the  study  of  the  sources  ship  structural 
damping  would  be  to  attempt  to  find  the  damping  effects  of  painting  and  lagging  on  the 
ship  structure. 
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